A quantum generalization of the thermal viscous friction law 
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Via 



On the basis of the equivalence of the energy balance de- 
scription at micro- and macro-level we propose a quantum 
generalization of the viscous friction law for a macroscopic 
Langevin equation describing thermal fluctuations without 
the zero point contribution. This equation recovers the clas- 
sical case in the limit of h — * 0. In particular it satisfies the 
quantum regression theorem and resolves several anomalies 
appearing in the quantum extension of the fluctuation dissi- 
pation theorem. 

PACS numbers: 03.65.Ca, 05.30.-d, 05.40.+j 



The classical Langevin equation for the damped har- 
monic oscillator with eigenfrequency oj s in the presence 
of friction and of a stochastic Langevin force fc(t) in the 
time domain writes: 
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dt 



dt 
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Here 7 is the viscous friction coefficient and the correla- 
tion function of the thermal Langevin force is given by 
< fc(0)fc(t) >= 2kTjS{t). Equation (1) forms a ba- 
sis for the fluctuation dissipation theorem (FDT) in its 
classical form: 
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(2) 



with S^x the two-side spectral density of the displace- 
ment fluctuations in the frequency range [— 00, 00] and 
a^(u) the classical susceptibility (i.e. the Fourier trans- 
form of the response function) that for the damped har- 
monic oscillator writes: 
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According to a recent comment of van Kampen 1 : the 
shortcoming of the above Langevin scheme is that it can- 
not be extended to quantum mechanics; attempts have 
been made to simple generalize Eq. (1) so as to turn 
it into an equation of operators, but that did not work. 2 
Indeed, one way to perform such a generalization is to 
introduced a relaxation already at the level of the Heisen- 
berg equations of motion for operators, the so called 
quantum Langevin equations. 3,4 However, as mentioned 



by van Kampen, such an approach is internally contra- 
dictory since formally it implies that a set of quantum- 
mechanical operators and their commutation rules which 
describe the physical system are relaxing. This would be 
a severe violation of the laws of quantum mechanics since 
commutation relations must be valid at all times. 4 

To overcome this contradiction, a second way of pro- 
ceeding is to introduce a relaxation for observables only. 
Accordingly, relaxation appears just in the macroscopic 
equations of motion for the observables, such as x — 
Tr{px}, that are obtained after averaging over the statis- 
tical operator /5. 4,5 This second approach is free from the 
above contradiction since the quantum corrections to the 
Langevin scheme are formulated only at the macrolevel 
on the basis of the quantum FDT (QFDT) given usually 
by the well known Callen and Welton 6 expression. In full 
analogy with the classical case, here the friction remains 
viscous and the modifications concern with the Langevin 
force spectrum only. However, some contradictions ap- 
pear in this scheme too. For example, when thermal 
fluctuations are considered this second approach enters 
in conflict with the quantum regression theorem (QRT) . 
Usually such a conflict is announced in the form "there 
is no quantum regression theorem" 7 ' 8 , despite of the fact 
that the QRT was derived by Lax 9 independently from 
the QFDT. 

The aim of this letter is to remove these contradictions 
of the second approach by introducing a quantum viscous 
force law which generalizes the relaxation of observables 
given by Eq. (1) to the quantum case. To this purpose: 
(i) the classical viscous friction law is replaced by the 
quantum analog as 
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and (ii) the classical Langevin force correlator by the 
quantum analog: 



</o(0)/Q(t) >=^Pe*p(-^ 
7i"P Jo 2 



)cos{ujt)duj (5) 



where 1/0= kT/h. 

The above quantum generalization of the Langevin 
scheme satisfies the QFDT in the form originally pro- 
posed by Nyquist 10 without the contribution of the zero- 
point energy of the thermal bath radiation: 
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Sg(cu) = 2hsgn{u}N T {u)Im{c%{u>)} 
with the quantum susceptibility a^(cj): 



(6) rewrite Eq. (8) as: 



ujI — u 2 — i^sinh(^) 



(7) 



Here Nt(uj) = [exp(@\u>\) — is the average number 
of thermal photons with frequency uj. 

Below, we shall demonstrate that the above extended 
approach is consistent with the principle of detailed en- 
ergy balance, while the conventional approach based on 
the classical viscous force and Callen-Welton expression 6 
of the QFDT violate this principle. 

Micro- and macro-level of detailed energy balance de- 
scription. We shall consider a sufficiently large isolated 
system, with Hamiltonian H = Hg + Hj- + V, which can 
be decomposed into two interacting subsystems. The first 
one is described by Hs and corresponds to the physical 
system under test. The second subsystem is described 
by Ht and represents the surrounding world. The in- 
teraction between these subsystems is described by the 
Hamiltonian V = —xf, where x and / are operators rep- 
resenting dynamical variables pertaining to the S- and 
T-subsystems, respectively. 

By using the standard procedure 4 to go from the den- 
sity matrix of the whole system p to the reduced density 
matrixes of subsystems ps = Trx{p} and px = Tr${p\ 
in the framework of the energy representation of pi — 
Pi(Hi) one obtains the following equation for the time 
evolution of the average energy < Hi > in the i-th sub- 
system (i = S,T): 



d £j d £r 

— < H T > = --77 < H s >= 
at at 
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From Eq. (11) we obtain the condition of macroscopic 
detailed energy balance (MaDEB) as: 



J xx {-U)jff{uj) = J xx (u)Jff(-U}) 



(12) 



which requires to be fulfilled for any value of the current 
frequency uj. We notice that the condition given by Eq. 
(12) is not the only form in which the MaDEB can be 
expressed. By using the expression for the imaginary 
part of the macroscopic susceptibilty 12,13 : 



Im{ai(uu)} = ^[JiiM ~ Jii(-u)] 



(13) 



which describes the dissipative part of the response of one 
subsystem to the action of the other one (for example, 
if i = x is the response then the action is given by / = 
Tr{prf} and vice versa) the MaDEB conditions can be 
rewritten in an equivalent form as: 

Y£(u,)Im{a f (u>)} = Y?{u>)Im{a x (u)} (14) 

where, in general, the weighted symmetric spectral den- 
sity Yf(cj) is given by: 



Y/» = (1 - /i)J«(-M) + nM\u>\) 



(15) 
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where p l m is the probability to find the i-th subsystem 
in the state with eigenenergy E l m , ui l mn = (E l m - ED/% 
the frequency of permitted transitions, x mn and Jmn the 
matrix elements of the x and / operators, respectively. 
When < Hi >= const, from Eq. (8) one directly obtains 
the condition of microscopic detailed energy balance 11 
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To formulate the conditions of energy balance at the 
macroscopic level of description, we replace the factor 
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using the definition of spectral density given by: 

Pn\^mn 



(10) 



which corresponds to the quantum correlation function 
Tr{pi(±T)i(0)} of the dynamical operators i — x,f, we 



and the weighting factor p can be treated as an arbitrary 
parameter. Note, that to derive Eq. (14) one need to 
preserve the antisymmetric property of the term in the 
square brackets of Eq. (11) under all the transformations. 
Since by definition Im{ai(uj)} is an antisymmetric func- 
tion of frequency, to keep the antisymmetry of the whole 
expression, l7*(w) must be treated as a symmetric func- 
tion, i.e. as a function of the absolute value of frequency. 

In the following we shall consider the formulation of 
MaDEB conditions by using the ratios of some charac- 
teristics of only one of the two subsystems. In such a 
representation, Eq. (12) takes the form: 



Jxx{u) Jff(u) 



p(u>) or Ju(-uj) = p(lo)Ju(lo) 

(16) 



where p(u>) is a factor common to both subsystems sat- 
isfying the relation p(—w) = p _1 (w). The MaDEB con- 
dition given by Eq. (14) can be represented in the anal- 
ogous forms: 



Yf(u,) = <f(u>)Im{ ai {u>)} 



(17) 



where 
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Here, g^(co) depends on p and is determined by the fre- 
quency dependence of p(u) at u > only. 

When the physical system interacts with the radiation 
bath described by a certain, not obligatory thermal, dis- 
tribution of photon numbers in the radiation modes, the 
bath spectral densities J//(±o;) can be represented as 14 : 



Jff(±\ 
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where G(u>) is the density of radiation modes in the bath 
and 7 2 (w) the electro-dipole matrix element square. Sub- 
stitution of Eq. (19) into the MaDEB condition given by 
Eq. (16) allows us to rewrite the latter in terms of the 
average number of photons in the radiation mode with 
frequency ui: 



(20) 



Accordingly the MaDEB condition given by Eq. (17) 
takes the form: 

Y£{u) = 2hsgn{uj}[N(uj) + p]Im{a x {uj)} (21) 

By using Eq. (21), let us give an interpretation to the 
energy balance at the macro-level. As follows from Eq. 
(14) the energy balance description at the macro-level is 
formulated only in the framework of the notion of a clas- 
sical absorption (i.e. the energy dissipation by one sub- 
system from another) which is described by Im{ai(uj)}. 
With respect to each interacting subsystem such an en- 
ergy transfer occurs in one direction, from outside to in- 
side. The opposite process, which returns the energy 
back, i.e. the energy emitted by the subsystem is absent. 
This representation is due to the definition of a>i(u)) as 
the response of the system to a classical external force 
which, by definition, does not take into account the pro- 
cess of spontaneous emission. Therefore, with respect to 
a single subsystem, one side of Eq. (14) must be treated 
as the absorbed (dissipated) power and the other side as 
the emitted (returned) power. Indeed, a direct substitu- 
tion of Eq. (19) into Eq. (14) gives us again Eq. (21). 
Here, the right-hand side, which contains Im{a x (Lo)}, 
describes the power gained from the radiation bath and 
dissipated in the system. The left-hand side, which is 
given identically by YJ?(uj), describes the power returned 
by the system to the bath. However, as follows from Eq. 
(21), at the MaDEB level the state of the bath radiation 
which is proportional to (N(u>) + p) and the spectrum 
of fluctuations Y£(uj) are formally uncertain since p can 
take arbitrary values. 

Equivalence between MiDEB and MaDEB. As follows 
from the above consideration, micro- and macro-levels 
of the energy balance description can not be treated as 
entirely equivalent. Indeed, the former is formulated in 



terms of statistical operators pi (see Eq. (9)), which de- 
scribe the interacting subsystem microstates only. By 
contrast, the latter is based on characteristics which are 
by definition invariant with respect to these microstates 
since all the MaDEB conditions are formulated for the 
quantities statistically averaged over fa. Moreover, the 
MaDEB conditions can include the extra arbitrary pa- 
rameter p which does not directly follow from MiDEB. 

It is easy to see, that the energy balance descriptions 
based on Eqs. (8) and (11) will be equivalent if the follow- 
ing relation is satisfied at all the frequencies of permitted 
transitions between the subsystems: 
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(22) 



At thermal equilibrium, when pi ~ exp(—Hi/kT), the 
common factor p(uj) — exp{— Tiuj/kT) becomes the uni- 
versal function of the temperature and frequency only, so 
that the MaDEB condition given by Eq. (21) at p = 
and 1/2 coincides with the Nyquist 10 and Callen-Welton 6 
versions of the QFDT, respectively. This fact allow us to 
give to the QFDT an alternative physical interpretation 
with respect to the conventional one as a relation which 
describes the MaDEB between two interacting physical 
systems under thermal equilibrium. 

The explicit presence of the uncertainty introduced by 
p ^ in the MaDEB given by Eqs. (17) and (21) leads to 
a disagreement between the micro- and the macro-level of 
description since, in essence, it claims that the statistical 
and photon-number operators, pr and n, respectively, are 
not sufficient to determine the average number of photons 
in the mode, i.e. N ^ Tr{prn}. In accordance with the 
MiDEB and MaDEB equivalence condition given by Eq. 
(22), p can always be neglected by renormalizing the av- 
erage number of photons in the bath N'(lu) = N(uj) + p. 
Such a renormalization will change only the value of ps, 
i.e. a microstate, by transforming YJf(uj) to J xx (— |w)|. 
In this case, only the photon part, HloN(uj), of the full en- 
ergy of the field is involved in the dissipation process de- 
scribed by Im{a x (u>)}, and the process of returning back 
the energy is accomplished by the spontaneous emission 
characterized by the spectral density J xx {— M). Thus, 
when p = 0, Eq. (21) describes the energy exchange 
between the two subsystems consistently with the micro- 
scopic picture. 

As follows from Eq. (22), in the frequency regions 
where p{w) = 1 the energy transitions are under satu- 
ration, that is p m — p n at u> = Lu mn . The saturated 
transitions can only appear jointly in both the interact- 
ing subsystems. In the opposite case, the MiDEB condi- 
tions are violated. By using Eq. (13) and Eq. (16), the 
imaginary part of the susceptibilities can be represented 
as: 

1 



Im{on{u)} 
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[1 - p(u>)]Ju(u>) 



(23) 



As follows from Eq. (23), at the frequencies of the sat- 
urated transitions the macroscopic response of both the 
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subsystems must not contain the dissipativc component, 
since Im{ai(cu)} = at p(u>) = 1. 

If one of the subsystems is the radiation bath, then, as 
follows from Eq. (20), the dissipativeless interaction oc- 
curs at all the frequencies where N(v) — ► oo. Note, that 
the state of the radiation in these modes is similar to the 
classical state with exactly determined phase 14 . Under 
thermal equilibium the frequencies of the saturated tran- 
sitions are determined by the condition exp{—[3uj) = 1 
and they correspond to the so-called Matsubara frequen- 
cies oj = Ofe = i2iTf3~ 1 k with k = 0, ±1, ±2, ... placed at 
k in the imaginary axis. The classical viscous friction 
law introduced usually to describe the macroscopic dissi- 
pation satisfies the MiDEB conditions of the dissipative- 
less interaction, however, on the real frequency axis only. 
Indeed on the real axis N(uj) — ► oo only at u> = where 
the viscous friction law always leads to Im{a x (oj)} = 0. 
However, this law fails by cxtendinding these MiDEB re- 
quirements to the whole complex frequency plane since 
it formally violates the energy balance at the Matsubara 
frequencies with k ^ where Im{a x (Cl] e )} 0- Such a 
violation is just the source of the QRT-QFDT conflict 7,8 
and of unremovable divergences of the correlation func- 
tions of fluctuations when T — > 0. 3 ' 5 

The Langevin scheme. In the framework of linear re- 
sponse theory, the spectral densities of fluctuations of 
the observable and the external force by which they are 
initiated are interrelated by the simple relation: 



(24) 



This allows us to rewrite MaDEB condition given by Eq. 
(17) in the form: 



Y?{u,) = -g»{w)Im{a- 1 (u)} 



(25) 



where Im{a~ 1 (u))} = — Im{ai(u))}/\cti(u))\ describes 
the relaxation law of fluctuations of the observable. For 
instance, Im{a~ 1 (u))} = —7a; in the case of the classical 
viscous friction. In the framework of the Langevin ap- 
proach the quantity Y^(u>) determines the spectral den- 
sity of the Langevin force which initiates fluctuations of 
the observable described by the spectral density YJf(u>). 

To match the Langevin scheme with MiDEB it is nec- 
essary: 

(i) To exclude all the extra contributions caused by the 
free parameter /1. It means that the spectrum of fluctu- 
ations of the observable must coincide with the sponta- 
neous emission spectrum of the system under test. This 
MiDEB requirement agrees with the opinion often meet 
in the literature 15,16 that zero-point fluctuations which 
corresponds to MaDEB condition with \x = 1/2 cannot 
be directly detected. 

(ii) To satisfy the dissipativeless interaction require- 
ment, the condition Im{a x (w)} = must hold in all the 
points of the complex frequency plane where p(ui) = 1. 

The former requirement is satisfied in a trivial way 
by assuming \i — 0. One of the possible ways to sat- 
isfy the latter one is to modify the viscous friction law. 



For this sake, by taking into account that by definition 
Im{a~ 1 (uj)} and Yj(ui) are odd and even functions of 
frequency, respectively, let us rewrite Eq. (25), which 
determines the specral density of the Langevin force, in 
the form: 



7 _£M 

2h-e > 



p Jm{a-V)} 

7 (e 2 — e 2 ) 
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In the classical limit fi —> 0, the first factor in the r.h.s. 
of Eq. (26) gives the classical spectral density of the 
Langevin force, 2kT"f, while the second factor is identi- 
cally equal to unity for the case of a classical viscous fric- 
tion, i.e. when Im{a~ (u>)} = —jlu. By assuming that 
the second factor keeps the same value also at h ^ 0, one 
obtains the quantum generalization of the thermal law of 
viscous friction in spectral representation as: 



Im{a- 1 (u,)} = l(e-^-e^) 



(27) 



which satisfies the second requirement of MiDEB in the 
whole complex frequency plane. The corresponding time 
representation is given by expression (4). In so doing, 
the first factor in the r.h.s of Eq. (25) gives the quantum 
analog of the spectral density of the thermal Langevin 
force which corresponds to the correlation function given 
by Eq. (5). 

In conclusion, we have proposed a thermal quantum 
viscous friction law at a macroscopic level. This law 
leads to a macroscopic quantum Langevin equation that 
does not include the zero point contribution in the spec- 
tral density of the fluctuating observable under thermal 
equilibrium conditions. The most relevant implications 
of this generalization are: (i) in the classical limit h — > 
the quantum Langevin equation recovers the classical one 
as it should; (ii) it resolves several anomalies appearing 
in the quantum extentions of the fluctuation dissipation 
and regression theorems, such as the well known QRT- 
QFDT conflict, and various divergencies originated by 
the Callen-Welton form of the QFDT; (iii) the predic- 
tion of an exponential like decay of the susceptibility at 
uj > kT/h. 
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